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nonlinear response with fluctuations far from thermal equilibrium and hold only in case of chiral 
transport. They can be used to test chiral or non-chiral character of the system. 
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According to the causality principle, past events in- 
fluence the future but the future has no effect on the 
past. This principle has no general counterpart in terms 
of the spatial separation of events: consequences of some 
events can be felt in every point after a sufficient wait 
time. A spatial version of the causality principle emerges 
in low-energy effective theories of some many-body sys- 
tems. The best known example is the integer quantum 
Hall effect (QHE): low-energy excitations are confined to 
the edges and can propagate only clockwise or counter- 
clockwise [l| . This can lead to a situation in which earlier 
events affect only those future events that occur "down- 
stream" . Similar chiral transport is possible in a number 
of other systems: some fractional quantum Hall liquids 
[lj], interfaces of topological insulators, superconductors 
and ferromagnets @, 0] , surface states in 3D QHE and so 
on. The simplest example comes from the statistical me- 
chanics models of traffic Q: chiral transport is possible 
on a network of one-way roads as long as no traffic jams 
form. 

In this paper we explore consequences of the extended 
causality principle in chiral systems. Causality is crucial 
for linear response theory. One of its celebrated results 
is the fluctuation-dissipation theorem (FDT). We show 
that a family of generalized FDTs holds in chiral sys- 
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FIG. 1: (color online) (a) r reservoirs are connected to sub- 
system S by chiral edge channels 0,1, ... ,r. The propagation 
direction along each edge channel is shown with an arrow. 
Each reservoir is at equilibrium with its own temperature and 
chemical potential. Panel (b) illustrates a quantum Hall bar. 
Transport occurs along chiral edges and in quantum point 
contacts. Dotted lines show tunneling between edges in the 
point contacts. 



terns. While the usual FDT applies in thermal equilib- 
rium only, our theorems are also valid in non-equilibrium 
steady states. 

The simplest relation 0-0] of such sort was derived for 
the exactly-solvable chiral Luttinger liquid model with a 
single impurity. We have recently found an FDT-type 
relation between the current noise and nonlinear con- 
ductance in a general chiral system in a non-equilibrium 
steady state in a three-terminal geometry Q. In this 
paper we prove a much more general result: we express 
nonlinear responses of the currents of various conserved 
quantities, such as the electric current and thermal cur- 
rent, in terms of the second and higher order cumulants 
of the statistical distributions of the currents in a non- 
equilibrium steady state in a multi-terminal system with 
an arbitrary number of terminals. The generalization is 
achieved due to a much simpler approach. The result for 
the chiral Luttinger liquid model follows from its tech- 
nically difficult exact solution Q. A more general result 
[H was obtained with a simpler but still rather subtle 
method, generalizing the equilibrium Kubo formalism. 
In this paper we use a completely different trick based 
on fluctuation relations [t| [Io[ . 
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The fluctuation theorem has been used in Ref. 
derive universal relations for nonlinear transport coeffi- 
cients in the absence of time-reversal symmetry (see also 



Refs. [12H19|V The results of Ref. [Ill hold irrespec- 
tive of chirality. Our results apply to chiral systems only 
and can thus be used to test transport chirality exper- 
imentally. This problem is of great interest 2Cj-|24[ for 
the QHE physics, in particular, at the filling factor 5/2. 
Indeed, the question of chirality of the edge transport is 
relevant for the current search for non-Abelian anyons 
at that filling factor [2^, [22j . Even in better understood 
QHE states, such as the Laughlin series, there is no com- 
plete theory of edge transport. In particular, the chiral 
Luttinger liquid model [l[ faces difficulties (for a review 
of recent experiments see Ref. |25[ ) . Our theorems provide 
a way to test its basic assumption of chiral edge transport 
[H [23] and touch upon the problem of nonlinear transport 
in QHE lUEzl. 



Non-equilibrium fluctuation-dissipation relations have 
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been derived in various classical systems [28H39J . In con- 
trast to Refs. 28-3^, we consider chiral transport and fo- 



cus on topological quantum systems in some of which chi 



ral transport has received experimental support [211 . 122 



We begin with a proof of a relation between the non- 
linear conductance and the current noise (second order 
current correlation function) and then generalize it for 
higher order correlation functions. 

The system (Fig. 1) consists of a central subsystem S 
and r reservoirs which we number with k = 1, . . . , r. We 
assume that the reservoirs are connected to S by chiral 
edge channels as shown in Fig. 1. Subsystem S may but 
does not have to be chiral (see an illustration in Fig. lb). 
There is no bulk transport outside system S and reser- 
voirs. Chiral edges and r are fully absorbed by reser- 
voirs 1 and r respectively. In the case of a QHE system, 
this means that the resistance of the reservoirs is much 
lower than the Hall resistance of the QHE bar. We do not 
make assumptions about the strength of the interaction 
between the remaining reservoirs 2, . . . , r — 1 and chiral 
edges. We assume that transport along edge is uncor- 
rected with what happens on the other chiral edges in 
Fig. 1, i. e., all forces are short-ranged. For a QHE bar 
this implies screening of the long-range Coulomb interac- 
tion by a gate whose distance from the 2D electron gas is 
shorter than the distance between edge and the other 
edges. Edge must also be far enough from the other 
edges to prevent transport into them through the bulk 
from edge |4l| ■ The absence of such leakage current 
can be verified by testing the quantization of the Hall con- 
ductance after all side reservoirs 2, . . . , r — 1 are discon- 
nected. Each reservoir is maintained at its temperature 
Tfc = 1/ /3k and chemical potential /i^ . Unless all temper- 
atures and all chemical potentials are the same, the sys- 
tem is in a nonequilibrium steady state. In order for the 
chiral description to hold, Tk and Hk must be lower than 
the QHE gap. Let the electric current Ii and heat cur- 
rent Ji flow into reservoir i. Our main result is that the 
zero-frequency cross-noises Su = lim w _j.o(A ( w Ki( — 
and = lim w ^ (^i( w )^i(~ w )) 5 are related to 

the nonlinear responses of Ii and Ji to the electrostatic 
potential V\ and temperature T±: 
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The Boltzmann constant ks is set to 1 throughout the 
paper. 

Our method is connected with the approach of Refs. 
(42I [43}. We consider the following protocol: Initially, 
subsystem S and the reservoirs are decoupled. An in- 
teraction V(t) that allows particle and energy exchange 
with the reservoirs is turned on at times < t < T. The 



interaction is turned off at t > T. At t < 0, reservoir i 
is at equilibrium with an inverse temperature /3j = 1/Xi 
and a chemical potential fii = qVi, where q is the charge 
of a charge carrier and V* the electric potential. We use 
only one set of chemical potentials and thus assume that 
only one carrier type is present. The initial state of fi- 
nite subsystem S is irrelevant. It is convenient to regroup 
S (and the rest of grey area in Fig. 1) with one of the 
reservoirs [Z^, for example, the r-th reservoir. The inter- 
action V(t) becomes a constant Vo when fully turned on 
during r < t < T — r. We assume that t <C T and 7~~ is 
much longer than the relaxation time so that the system 
remains in a steady state during most of the time interval 

r. 

The above process is called a forward process in the 
formalism of fluctuation relations. We also need to study 
a backward process which can be described as a forward 
process in the time-reversed twin system with the op- 
posite chirality (i.e., the directions of all arrows must 
be reversed in Fig. 1; see Supplementary material for a 
detailed discussion). We assume that the initial temper- 
atures and chemical potentials of all reservoirs are the 
same in the forward and backward processes. 

Consider now the changes ANi = Ni(t = T) — N% (t = 
0) in the particle number JVj in each reservoir and the 
changes AEi = Ei(t = T) — Ei(t = 0) in the energy Ei 
in each reservoir. The total energy and particle number 
are conserved, so that J^i AEi = J2i ANi = 0. We in- 
troduce the joint probability distribution of the particle 
number and energy changes P[AE, AN;±], where the 
vector AE = {AEi}, AN = {ANi}, and the "+" or "-" 
sign in the last argument of P refers to the forward or 
backward process respectively. According to the fluctua- 
tion relation (a derivation can be found in Ref. [42[ and 
Supplementary material) 



P[AE, AN;+] 
P[-AE, -AN; -] 



n 



(3) 



Given a distribution P[AE, AN; v\ with v = ±, we are 
able to calculate the heat currents J" = limj-^oo^AEi — 
liiANi)) v /T and their correlation functions, as well as 
the electric currents I" = \im T ^ fOC (qAN i ) u /T and their 
correlation functions. The triangular brackets mean tak- 
ing average with respect to P[AE, AN;i/]. The limit 
T — > 00 is taken so that the steady-state quantities are 
obtained. It is convenient to define a cumulant generat- 
ing function 

Q(x, y, f3, n; v) = Yirn^ 1 In j / J] {dAE.dAN,) 

^ i<r 

x e -£i^ AS *-£^ A ^p[AE,AN;/3,/z;z;]j, (4) 

where the vectors x = {x.;}, y = {yi}, [3 = {A}, and 
fj, = {fJ-i}. The /x-dependence of Q cannot be reduced 
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to a dependence on the differences Hi — Hi even in the 
QHE context because of the screening gate. Since P is 
normalized, we have (3(0, 0, (3, fi; v) = 0. The currents 
and their correlation functions can be obtained by taking 
derivatives of Q(x, y,/3, /i; v) over or j/, and then set- 
ting x = y = 0. The n-th order heat current correlation 
functions 



C 



ij—k 



\y<lr /'<<>., <h HjBy 

■[d Xk - (i k dy k ]Q(0,0,f3,fi;v), (5) 



and the n-th order electric current correlation functions 



C 



(»)»" 

ij—k 



NP 9 A-^Q(o,o,A/x;4 (6) 



At n = 1, cM 1 " 1 ' and Cy~ hV are just the heat current 
J" and the electric current I" . The second-order correla- 
tion functions are the low-frequency noises S^- v = Cj^- 
and S"j = [this definition differs by a factor of 2 



from Ref. [8|. Higher-order correlation functions have 
also been studied experimentally (3]. The fluctuation 
relation (1A4[) leads to a symmetry of the generating func- 
tion 

Q(x,y,/3,M;+) = Q(/3-x,£-y,/3,M;-) > (7) 

with | = {-AmO- 

As discussed above, edge (Fig. f) is independent 
of the remaining edges. In other words, our protocol re- 
sults in two statistically independent transport processes: 
charge and energy transfer along the lower edge and 
along the remaining edges in the upper part of the sys- 
tem. This means that the distribution function can be 
rewritten as 

F[AE,AN ;i /]= / dAE[dAN[ Pi[AE[, AN[; v]P 2 [AEi 

-AE[, AJVi - AN[ , AE 2 , . . . A£ r _i, AN U . . . , AJV r _ i; u], 

(8) 

where Pi [AE[ , AJV{ ; v\ is the probability to trans- 
port AE[ energy and AN[ particles into reservoir 
I along the lower edge (edge in Fig. I), and 
P 2 [AE'l,AN[',{AE s },{AN a };v] is the probability to 
transport AE" energy and AJVf particles along edge 1 
and change the energy and particle numbers in reservoirs 
2, . . . , r - I by (AE^AJV,), I < i < r [Remember that 
AE r and AN r are not independent variables because of 
conservation laws]. Negative AE[, AE'{, AN[ or AN" 
mean the energy and/or particle loss by reservoir 1. We 
now discuss the dependences of P\ and P 2 on the temper- 
atures Pi and chemical potentials Hi- In the setup with 
the "+" chirality, energy and particles on the lower edge 
flow out of reservoir r and into reservoir I . Due to the ex- 
tended causality principle, the distribution of AE[, AN[ 
depends only on p r and fi r . Meanwhile, reservoir r re- 
ceives energy and particles from the upper part of the 



system but does not provide any feedback, so the trans- 
port in the upper part does not depend on p r and /i r . 
Hence, Pi only depends on p r and Hr while P 2 does not 
depend on P r and fi r . In the setup with the opposite "— " 
chirality, the same argument shows that P\ only depends 
on Pi and Hi while P 2 does not depend on Pi and Hi- 

In terms of the cumulant generating function, Eq. (8) 
means that Q(x,y,(3, fi; v) = Qi + Q 2 splits into two 
terms Qi and Q 2l corresponding to Pi and P 2 respec- 
tively: 



Qi(x,y,j9,/i;i/)= lim ^ In ( / dAE[dAN[ 
r^oo 7 l J 

, e -(xi —x r )AE' 1 — (y 1 -y^ANi Pj [ A _g/ ^ AJV j .p^-jA. (Q) 



Ki<r 



Q 2 (x, y, 0, (i; v) = lirn^ f^lj dAE'(dAN'{ 
dAE t dANiexp(-xiAE" - x i AE l - x r 5E r ) 

Ki<r 

x exp(-yiAN{' - ^ yiAN-y r SN r ) 

Ki<r 

xP 2 [AE';, AN[', {AE r>l>1 }, {AN r>i>1 }; f3, fi; v\ Jq.0) 

where we used the relations = S A Ni = and 

denned 6E r = -AE'{ - ^ 1<1<r AE„ SN r = —AN" - 
^2i <i<r AN. The chirality-induced causality means that 
Qi{v = +) depends only on p r and Hr while Q 2 {v = +) 
does not depend on p r and Hr- Qi(y = — ) depends only 
on Pi and Hi and Q 2 (v = — ) does not depend on Pi and 
Hi- 

We are now ready to prove the steady-state FDT for 
chiral systems. Let us start with the particle transport. 
We apply the differential operator £>i = D Vi — TiD^ H to 
both sides of Eq. ([7]). D Vi and D Mi stay for full derivatives 
over the respective variables. We obtain 

{d yi -T i d lli )Q(^,y,l3, l i;+) = -Tid^Qip - x,£ - y,/3, /x; 

(11) 

We emphasize that the partial derivative with respect 
to Hi on the right hand side is taken at a fixed The 
expression £j = —Pi Hi should be substituted after the dif- 
ferentiation. This reflects the difference of the operators 
D and d. We now apply the differential operator Dj on 
the two sides of Eqs. (|TT|) and set x = y = at the end. 
In terms of the correlation functions (9), one finds 



di+ bit 

°dV n + l dV, ~ « 



(12) 

In order to derive Eq. (fT2")) we use the iden- 
tity B^Bfj,, Q(0,0, (3, /x;+) = which follows from 
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Q(0,0,P,/i;v) = 0. Note that the last term in Eq. (QH) 
is defined for a system with the "— " chirality, whereas 
all other terms refer to the "+" chirality. 

In a chiral system the last term in Eq. (fT2|) is zero at 
i = 1 < j < r. This can be seen by writing the partial 
derivative ,<2(x = (3, y = £, /3, /x; — ) as the sum of 

the derivatives of and Q2- In the system with the 
"— " chirality, Qi depends only on fj,i while Qi does not 
depend on fix before we make the substitution x — > (3 
and y — > Thus, the partial derivatives d^d^Qi and 
d f i 1 d flj Q2 are both zero. Moreover, if 1 < j < r then the 
first term of Eq. (fT2|) dl± /d[ij is also zero since I± only 
depends on /3i jr and fj,\ tr in the system with the "+" chi- 
rality. Hence in systems with the "+" chirality, Eq. (fl2|) 
simplifies to Eq. (p}. This is our main result: The cross 
noise between the currents in reservoirs 1 and j is con- 
nected to the response of the current in reservoir j to the 
voltage in reservoir 1 regardless of the non-equilibrium 
nature of the system. If j = r then the term dI\/dV r 
is not zero. Instead, it equals the conductance G in the 
two-terminal setup with only two reservoirs 1 and r. This 
is the case since dI\/V r is solely determined by edge 
and hence does not change if all reservoirs 2, 3, . . . , r — 1 
are removed. Thus, at j = r there is an additional con- 
tribution GT r on the left hand side of Eq. (QJ. In a QHE 
system with the filling factor v, G — ve 2 /h. 

Similar results can be obtained for heat currents and 
noises after one applies the differential operators Dj = 
D Xj - LijD Vj + Dfy and D\ on both sides of Eq. 0. For 
the systems with the "+" chirality and j ^ 1, r, we derive 
Eq. (J3J) with this trick. For j = r, there is an additional 
term dJ\/d(i r = — (T r ) 2 dJi/dT r on the right hand side of 
Eq. @. In QHE, the physical meaning of the additional 
contribution d J\ /dT r is the thermal Hall conductance of 
a two-terminal Hall bar, equal to Kir 2 T r /3h, where k is 
universal. 

We stress that all above response functions are non- 
linear responses in non-equilibrium steady states. The 
FDTs ((T|) and @ hold for all chiral systems, as long as 
Ti and Vi are smaller than the energy gap in the bulk. For 
non-chiral systems, such as QHE systems with charged 
or neutral counter-propagating modes at the edges, the 
last term in Eq. (|12|) is nonzero in non-equilibrium states 
and the theorem does not apply. One microscopic mech- 
anism of FDT breaking in non-chiral systems involves en- 
ergy transport by 'upstream' modes from region S (Fig. 
1) to reservoir 1 along edge 1. Local heating at the hot 
spot, where edge 1 enters reservoir 1, affects the noise 
of the current, emitted from reservoir 1, as well as the 
cross- noises Su- 

Eqs. (fl"l2j) differ from the related results of Refs 0, B 
In the case of Ref. |7| this reflects a different geometry. 
Ref. [H considers a version of Fig. 1 with r = 3 and 
T r = Ti (Fig. lb). We verify below that the result of 
Ref. B can be derived from Eq. (|T|) with r = 3. Charge 



conservation implies that for low-frequency components 
of the electric currents I3 = —Ii — 1%. Hence, S33 = 
Su + S22 + 25*12, where Su denotes the auto-correlation 
noise and S12 is the cross-noise (P). Since the edges, 
connected to reservoir 1, are always in equilibrium with 
the temperature T\ = T r , we find Su = 2GT\ from the 
equilibrium Nyquist formula [remember a missing factor 
of 2 in our definition of the noise] . Then from the chiral- 
system FDT (JTJ) we get 

S , 33=S , 22-2T 1 ^- + 2GT 1> (13) 
aVi 

in agreement with Ref. 0. Certainly, the results of the 
present work go well beyond Refs. 0, H we cover multi- 
terminal geometries, thermal currents and a family of 
FDT's for higher-order correlation functions which we 
derive below. 

The strategy of their derivation is the same as above. 
We apply £>i on both sides of Eq. (J7J) m times and obtain 

b i --b j b h Q{o,o,i3, l x-+) = 

{-l^TiTj ■ ■ ■ T k d^ ■ ■ ■ d^ h Qip, t (3, /*; -). (14) 

We set m > 2, k = 1 and j 7^ l,r. The right hand side 
vanishes due to chirality. The left hand side expresses 
as a combination of derivatives of correlation functions. 
It simplifies dramatically after we combine Eqs. (fT4")l for 
all m < n. Specifically, one can prove by induction the 
following FDT: 

dc (n - 1] ' + dc {n -p' + d 2 c (n - 2h+ 

Ci-il Q V . Q Vl 

^ (15) 

where a hat above an index means that the particular 
index is absent in the index set (but does not mean that 
its value is absent in the set since we do not exclude a 
situation with two or more identical indexes). To derive 
relations between heat current correlation functions and 
their response functions, we apply the differential oper- 
ator &1 multiple times on both sides of Eq. (J7]). We 
obtain 

gQh(n-l) Q£,h(n-1) Q 2 (jH n - 2 ) 
r< h ( n ) j_(rn \2 , (rp »-'-ji , (rp rp \2 _ 

(T ' } &Tj +(Tl) m HTjTl) dT.dT, 

(16) 

Eqs. (|l|2p are special cases of (fT5|) and (fTBf . 

In conclusion, we established a family of fluctuation- 
dissipation theorems for charge and heat transport in 
chiral systems in nonequilibrium steady states. The ar- 
gument combines the formalism of fluctuation relations 
with the extended causality from chirality. Our results do 
not hold in nonchiral systems away from equilibrium and 
can thus be used to probe the chiral character of charge 
and energy transport. This question is of significant cur- 
rent interest [20M24| in the QHE physics and is relevant 
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in other fields, such as transport in heterostructures of 
topological insulators. 
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Supplementary material for Chirality, causality, and fluctuation-dissipation theorems in non-equilibrium 

steady states: Derivation of the fluctuation relation 

To make the paper self-contained we include a brief derivation of the fluctuation relation (|A4j) [Eq. (3) of the main 
text]. Additional details can be found in Refs. [42,43] and references therein. 

The main text defines forward and backward processes. We want to connect the joint distribution functions 
P[AE, AN; ±] of the energy and particle number changes in the reservoirs in the forward and backward processes. 
We first derive an expression for the distribution function in the forward process. 

We want to find the statistical distribution of the changes ANi = Ni(t = T) — Ni(t = 0) in the particle number 
Ni in each reservoir. Let Hi and Mi be the Hamiltonian and particle number operators of the i-th reservoir (H r 
includes system S). The particle numbers conserve in the absence of V(t), i.e., [Hi, Mi] = 0. Thus, the initial density 
matrix factorizes into a product of Gibbs distributions in each reservoir, p n = \\ { e -Pi[ E in-mN in \ /Zq, where Zq is the 
initial partition function and the index n labels the quantum state \ip n ) with the reservoir energies B m and particle 
numbers N. in . An initial joint quantum measurement of Hi and Mi is performed at t = 0, so that the quantum state 
of the system collapses to a common eigenstate |-0n) with the probability p n . Then the state \ip n ) evolves according 
to the evolution operator U(t; +) determined by the Hamiltonian H(t; +) = Ei % + ^W- At £ = T, a second joint 
measurement is taken, leading to the collapse of the system to the state \ip m ) with the reservoir energies Ei m and 
particle numbers N im . The probability to observe such process is P[m, n] = \(ip m \U(T; +)\ip n )\ 2 Pn- Hence, the joint 
distribution function of the energy and particle changes AEi^ mn = Ei m — Ei n and ANi^ nn = Ni m — Ni n is 

P[AE, AN;+] =^2Y[5(AEi - AEi, mn )S(ANi - ANi, mn )\(^ m \U(T; (Al) 

rnn i 

where the vector AE = {Ai^} and AN = {AiV;}. The total energy and particle number are conserved, so that 
Ei AE i = Ei AN i = °- 

We also need to study the backward process which can be described as the forward process in the time-reversed 
twin system with the opposite chirality. The time evolution operator U(t; — ) of the twin system is determined by 
the Hamiltonian T-L(t;—) = OH(T — t;+)Q~ 1 , where 9 is the time-reversal operator. The i-th reservoir has the 
Hamiltonian OT^G -1 . Clearly, 9|V>m) is a common eigenstate of QTLiO^ 1 and Mi = QMiQ^ 1 with the eigenvalues 
Ei m and Ni m . We assume that the initial temperatures and chemical potentials of all reservoirs are the same in 
the forward and backward processes. Note that such initial temperatures and chemical potentials are not necessarily 
the same as the final thermodynamic parameters in the forward process for large but finite reservoirs. Performing 
two quantum measurements at t = and t = T in the beginning and end of the backward process, one finds the 
distribution of the energy and particle number changes, similar to Eq. (|Aip : 

P[AE, AN; ] = II S ( AEi - AE ^m) 

mn i 

x S(AN t - AA^ nm )|(Vvi|0 _1 C/(T; -)e\^ m )\ 2 Pm , (A2) 

where the initial density matrix p m = Yii e~ ,3i ^ Eim ~^ iNim \ /Z^ . It follows from the antiunitarity of G that Zq = Zq . 
The evolution operators have an important property [41] 

Q- 1 U{T;-)Q = U\T;+). (A3) 

Combining Eq. (|A1[) . (|A2|) and (|A3[) . we recover the fluctuation relation 



P[AE,AN;+] =n » (ABr(liMi) 
P[-AE,-AN;-] 11 



(A4) 



